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04 . We investigate superconducting proximity effect in clean ferromagnetic layers with rough bound- 

aries. The subgap density of states is formed by Andreev bound states at energies which depend 
on trajectory length and the ferromagnetic exchange field. At energies above the gap, the spectrum 
is governed by resonant scattering states. The resulting density of states, measurable by tunneling 
spectroscopy, exhibits a rich structure, which allows to connect the theoretical parameters from 
experiments. 
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Investigating of sujperconducting proximity effect in normal systems has a long history back to the experiments of 
McMillan and RowellEJ. Their tunneling spectroscopy measurements in normal metals connected to a superconductor 
revealed strong modifications of the density of states (DOS) caused by the induced superconducting correlations. 
These results were understood in the tunneling model of McMillana. He first noted that the changes in the DOS of 
the normal metal occur on a scale £Th, which in his model is identified with inverse escape time of a quasi-particlc 
in the normal metal. Naively, one would have expected A, the superconducting gap parameter, to play the dominant 
role, which turned out to be not the case. Furthermore, if -&rh <C A, the energy gap in the superconductor plays no 
role at all and the DOS in the normal metal has a gap of the order of -&rh- We can understand these observations 
by noting that in clean normal metal films the electronic properties are determined by so called Andreev bound 
t— ( \ statesa. These are bound electron-hole pairs living on trajectories, which start and end at the superconductor. The 
J> ■ characteristic energy scale in this case is vf/(1, again the inverse escape time. In recent years experiments became 
possible, in which the density of states was resolved locally on a sub-/iin scale. For example, the dependence_of the 
tunneling DOS on the distance from the superconductor in normal metals has been measured by Gucron et al.u using 
additional tunnel junctions. These results have been successfully explained within the quasi-classical theory in the 
diffusive limit B Other experiments made use of low temperature scanning tunneling microscopes to resolve spatially 
the DOS of small droplets of normaLrnetal on the surface of a superconductor.! Nowadays these types of experiments 
^ , are becoming a standard technique J3~u 

The question of the proximity effect in the presence of a spin splitting is currently heavily investigated. In particular 

S the influence of a superconductor on transport properties of a ferromagnet is under debatell3~113,in which case the 
proximity effect is negligible. It is, however, natural to address the question of the influence of an exchange field on the 
proximity density of states. In fact, this question was already addressed experimentally a while ago by Gallagher et 
al. l3. They observed a spin splitting of the DOS in thin normal layers in a parallel magnetic field. New experimental 
developments]!.] exploring the proximity effect on a nanometer scale demand new theoretical models, beyond the 
£j , simple tunneling model of McMillan. In the present paper we investigate a new model, suitable for these experiments. 

The motivation stems not only from fundamental question of the coexistence of ferromagnetism and supercon- 
ductivity, but also because interesting applications of ferromagnet-superconductor (FS) hybrid structures have been 
proposed. We only mention here the potential use of SFS-contacts in the construction of quantum computers. . SFS - 
junctions are candidates for all-electronic 7r-junction, which are needed in some proposals for solid-state qubits.F 5 !! 16 ! 

We also note, that many surprising effects in normal metal-superconductor heterostructures have been found, both 
experimentally and theoretically. The question, how these effects are modified by the presence of ferromagnetism, 
is of high interest. For example, making use of mechanically controlled break junctions, the conduction channel 
content of single atoms has been determined!!!!. These were shown to depend on the chemistry of the atom onlyO. It 
was however crucial to have superconducting contacts in these experiments. Therefore, so far only superconducting 
elements or normal metals, with the help of the proximity effect, have been investigated.L3 It is reasonable to expect 
that experiments with ferromagnetic materials will be performed_in the near future. 

Long time ago Larkin and OvchinnikovE3 and Fulde and FerreBJ independently predicted, that in the presence of an 
exchange field h (for instance in magnetic superconductors) a specific superconducting state can be formed, in which 



Cooper pairs have a non-zero momentum due to spin splitting. The momentum of Cooper pairs is given by 2h/v-p, 
which, in the singlet state, results from the difference between up and down spin band Fermi momenta. The resulting 
LOFF state is qualitatively different from the homogeneous zero-momentum state. Due to a varying superconducting 
phase, the order parameter has an oscillatory spatial modulation containing nodes, in which the phase changes by 
7r. The LOFF state has never been observed in bulk superconductor, but there are recently evidences for detecting 
an induced LOFF state in heterostructures of ferromagnets and-, superconductors. Many works have investigated 
the thermodynamics properties of FS-multilayers. Radovic et alE3 have predicted oscillation of the superconducting 
critical temperature T c as a function of the. thickness of the attached F-layer. The experimental evidence for these 
T c oscillation^ is not, however, conclusive .E3 The reason for this may, for example, result from a bad quality of the 
FS-interfaceEJ ■_. 

The most recent experiments have concentrated on other properties of FS-layers. Ryasanov et alrB measured 
the temperature dependence of the critical current in SFS Josephson junctions with thin F-layer and have found a 
non-monotonic temperature dependence. This behavior can be understood in terms of a 7r phase shift due to the 
exchange field, which occurs for certain values of the thickness of the F-layer, as-first was predicted by Bulaevskii et 
al.t3. An indirect proof of the 7r phase shift has been made by Kontos et a/.JiJ, who studied the density of states 
in thin ferromagnetic films contacted by a superconductor. They observed an oscillatory behavior of the induced 
superconducting correlation for layers of different thickness, which was attributed to influence of the exchange field. In 
Ref.E3, we have shown that these experimental findings could be explained by a model of a ballistic ferromagnetic layer 
with rough boundaries. The best agreement was obtained in the limit of large h/A and small interface transparency 
T. 

In the present paper, we study the proximity DOS in a clean ferromagnetic layer on top of a superconductor in 
the full parameter range. Within the ballistic quasi-classical formalism we obtain, that the DOS for energies below 
the superconducting gap A is completely specified by the length distribution of the classical trajectories inside the 
ferromagnet (Sec. ph. The DOS for energies above the gap is also expressed in terms of the length distribution of 
the trajectories. The length distribution depends on the geometrical properties of the attached ferromagnet and the 



connecting boundaries. In Sec. Ill we specify the classical length distribution for our particular case of the F-film 
geometry depicted in Fig. fil. We assume, that the boundaries of the F-film are disordered, leading to complete diffusive 
reflection of the quasi-particles from these boundaries. We also take into account band mismatch and disorder at 
the FS-interface, which leads to an enhanced backscattering from this interface. For simplicity, we assume a single 
value of the FS-interface transparency T. With the calculated distribution, the DOS at all energies is obtained as a 
function of the superconducting gap A, the exchange field h, the thickness of the F-layer d and the transparency T. 
We analyze the DOS for different regimes of h/A. It shows the interplay between ferromagnetism and superconduc- 
tivity depending qualitatively on the thickness d. For example, a weak exchange field leads to a spin-splitting of the 
DOS, which results in a distinctive low energy peak in the total DOS. In addition there is an overall suppression of 
the superconducting features of the DOS with increasing h. At higher exchange fields, the DOS shows as a signature 
of the exchange splitting an-oscillatory behavior as a function of the layer thickness. This oscillation of DOS was 
observed in the experimentalj. Our findings are summarized in the following list: 



h = (Sec. |VA|): Andreev levels are governed by the distribution of trajectory lengths, which only depends on 
the geometric properties of the sample. At small energies E -C vp/d the DOS is strongly suppressed, originating 
from the exponential suppression of long trajectories. The DOS at larger energies reflects the length distribution. 
In our model it display a multiple peak structure, resulting from multiple reflections at the SN-intcrface. The 
resulting minigap correspond to the gap found in a calculation including impurity scattering.Ea 



• 



h/A < 1 (Sec. IVB): A small exchange field 'splits' the DOS for spin up and down quasiparticles, i. e. the total 
DOS is more or less a superposition of 'normal' DOS's at energies E ± h. Accordingly, the former minigap in 
the DOS is destroyed. It only remains a dip in the DOS shifted to finite energies. The density of states at the 
Fermi level approaches the normal state values in an oscillatory way, i. e. overshooting the normal DOS for 
certain values of h. 



h/A > 1 (Sec. IV C ): The superconducting features of the DOS are stronger suppressed. The former peaks 
at ±A are inverted into dips for thicker layers. Above the gap peaks at E = ±/i appear as the signature of 
resonant transmission through the ferromagnetic film. For thin layers features at ±ft, are absent and the DOS 
approaches a BCS-form. 



h/A 3> 1 (Sec. 1VD): For layers with d>vp/h the DOS exhibits coherent oscillation, i. e. the form of the DOS 
difference from the normal state value becomes independent of d. The amplitude and sign, however, depends 
on the thickness. Only for very low thicknesses d <C Tvp/h the DOS approaches the BCS form. 



In Sec. |lVE we condense our results into a map of the proximity DOS. Finally we present some conclusions in 
Sec. 0. 



II. MODEL AND BASIC EQUATIONS 

The system we study is sketched in Fig. (jl. A ferromagnetic layer (F) of thickness d is connected to a superconducting 
bank (S) on one side and bound on the other side by an insulator or vacuum. F is characterized by an exchange 
splitting, which we take into account as mean field h in the Hamiltonian. The thickness d is larger than the Fermi 
wave length Af and smaller than the elastic mean free path £i mp , which allows for a quasi-classical descriptions in 
the clean limit. We apply the Eilenberger equation in the clean limit 



-zv F V<? CT (w,v F ,r) = \{E + ah(r))f 3 - iT 2 A(r),g (T {uj, v F ,r)] . 
The matrix Green's function for spin a has the form 
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It depends on energy E, the direction of the Fermi velocity v F and the coordinate r. Here fj denote the Pauli matrices, 
A(r) is the superconducting pair potential (taken as real), and a (= ±1) labels the electron spin. The matrix Green's 
functions obey the normalization condition cfe = 1. Inside then F-layer h is constant and A = 0. We neglect a 
depression of the pair potential close to the FS interface, thus A(r) = const, inside the superconductor, which applies 
in the case of a bad contact between the ferromagnet and the superconductor. Strictly speaking, we would have to 
include an elastic collision term in (pi), even in the limit £np P ^> d. However, changes in the spectrum due to this term 



are limited to small energies ^Svp/£i mp <C min(vp/d, /i)E3, which are negligible in all cases we study, except for the 
case h = 0. Disorder in the superconductor can be neglected in the limit of small interface transmission, which we 
mostly assume. 

We have to solve Eqs. (|l]) along each classical trajectory with length I in F, that comes from the superconductor 
and ends there. As boundary conditions the solutions approach the bulk values of g a at the beginning and, the end 
of a trajectory deep inside the superconductor. These are given by g CT (bulk) = (— iEfa + Ati)/VA 2 — i? 2 .E!l It turns 
out that on a trajectory inside F the normal Green function g a is constant. It depends only on the length of that 
trajectory I and is given by 



g a = tanh [(— iE — iah)l/v-p + arcsin (—iE/A)} 
To find the density of states per trajectory, we have to calculate 



N(E,l) = -£ J2 Reg CT (£ + zO,v F ,r). 



(3) 



(4) 



a=±l 




FIG. 1. Schematic drawing of our model of a ferromagnetic film (F) topping a superconductor (S). Typical classical trajec- 
tories are also indicated. We distinguish two processes. At energies E below the superconducting gap A, quasiparticles are 
confined to the film by Andreev reflection (indicated by the white circles). An examples for this process is given by the left 
trajectory. For E > A Andreev reflection is incomplete and the quasiparticles states in F are formed by scattering resonances, 
symbolized in the right process. 



As the result we obtain for energies below the gap (\E\ < A) 
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where 



(nn + arccos(£'/A)) . (6) 



E + ah 



Above the gap (\E\ > A) we find 



w = fE E u^^,^ M l^L^,^ - ^ 



(7 — ±1 71: 



-^ [(£ + cr/ l )Z/w F - nvr] 2 + (acosh|£/A|r 



Here Ao is the density of states at the Fermi level in the normal state. Eqs. (0) express the fact that the density of 
states below A is a sum of <5-peaks resulting from Andreev bound states of electrons of E > (positive n's) and holes 
of E < (negative n's). The energies also follow from the quasi-classical quantization condition I = l n . 

The total DOS is obtained by averaging the expressions (0) and (0) over all classical trajectories. Denoting the 
trajectory length distribution by p(l) and using (H), we find for the sub-gap DOS 

N(E) = / ' dlp(l)N(E,l) 



a— ±1 ' ' n— — oo 

This formula is a general result for the sub-gap density of states of a quasi-ballistic metal connected to a super- 
conductor. It is completely specified by the length distribution of classical trajectories, which depends only on the 
geometrical properties of the attached ferromagnet and the surrounding boundaries. 
Averaging expression (Q) over p(l) the total DOS for energies above the gap is 

^)4E E IW) UW „„, aC0Sh| f /A ' ^ for |^| > A, (9) 



cr— ±1 n— — oo 



[{E + ah)l/v F - H 2 + (acosh|S/A|) 



The absence of discrete bound states reflects the fact, that the Andreev reflection at energies above the gap is 
incomplete. Therefore, the quasi particle states in the ferromagnet are determined by 'scattering resonances' of 
quasiparticles incident from the superconductor. 

III. DISTRIBUTION OF THE TRAJECTORY LENGTH 

Now we specify the length distribution for our particular case. We model the F-layer by a weakly disordered thin 
film bounded by a rough surface to the insulator and a rough FS-interface of average transparency T. A typical 
classical trajectory is depicted in Fig \u. An electron coming from the bulk of S enters into the F-layer and after 
several reflection from the insulator and the FS-interface returns to the S-bank, where it is Andreev reflected as a 
hole and transverses the trajectory in the opposite direction. Thus, the building block of a trajectory is the segment 
between two successive reflections from S. The number of blocks which form the total trajectory depends on the 
transparency of the interface, i. e., it is roughly ~ 1/T. 

As first step, we consider the length distribution in the case of a perfectly transparent FS-interfaces, in which the 
length distribution is that of one elementary block. Due to the roughness of the insulator and the FS-interface the 
quasi-particles undergo diffusive reflection from these boundaries. Incident and reflected directions are completely 
uncorrclated. Then, assuming an uniform distribution for the directions of vp, we obtain for the length distribution 
of one elementary block (corresponding to the case of T = 1) 

po(0 = / d(cosft) / d(cos9 r )S(l - -?— - -A-), (10) 

J J cos0i cos6 r 



where 9i{9 r ) denotes angle of the incident (reflected) direction with respect to the normal to surface of the insulator. 
To take into account the weak bulk disorder, we include a factor exp— l/£[ mp . This serves mainly to yield a finite 
average length of the distribution (Jin). In a purely ballistic layer with po(l) given by (|10|), the average length would 
logarithmically divergent. Taking this into account we obtain 



PO(0 = 7^2 
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(11) 



where C — E%{d/li mv ) (^2(2) = J\ dxexp (—zx)/x 2 is the exponential integral of order 2). 

In the second step, we connect the elementary building blocks, if the FS interface has a transparency T < 1. In 
determining the length distribution we assume that an particle either goes through the interface or is fully reflected. 
Only the number of these reflection depends on T. We do not take into account quantum mechanical interference for 
a single reflection at the FS-interface. Taking this into account will lead essentially the same results as our approach. 
By an expansion in the reflectivity R = 1 — T for the distribution p(l) we can write 



00 p n 

p(l) = Tj2R n dl ..dl nP o(lo)-Po(ln)S[l-Y / l i], 

n=Q J i=0 



(12) 



where the nth term in the expansion is the contribution of the trajectories on which quasi-particles after n times 
reflections from FS-interface leave the F-layer. It is easy to see from Eq. (^2|), that p(l) obeys the integral equation 



p(l) = Tp (l) + R / dl'p (l')p{l - I') 



which is readily solved by a Fourier transformation: 

p(l) = f^e ikl P(k). 

Replacing (|l|) in Eq. @ we find 

TPo(k) 



P(k) 



l-RP {k) 



(13) 



(14) 



(15) 



where Po(k) — E%(ikd + d/£[ mp )/C is the Fourier transform oi p (l). 

The distribution p(l) determines the relevant length scale associated with the geometrical size of the system corre- 
sponding to the typical distances quasi-particles travel inside F. We have plotted p{l) for different Ts and d/t lmv = 0.1 
in Fig. g, For small T, it has a characteristic double peak structure close to the shortest trajectories I ~ 2d, resulting 
from trajectories reflected once and twice from the insulator. At large I the distribution decays exponentially as 
exp(— I /I), where I ~ 2dln(d//i mp )/T is the mean trajectory length. We therefore have two characteristic lengths 
of the distribution, the smallest possible trajectory length 2d and the average length I. The former determines the 
energy of the first Andreev level and the latter the possible longest length of the trajectories. For T ~ 1 these two 
length scale are of the same order, leaving the thickness as the only relevant length scale. In this case p(l) has only 
one peak close to 2d (see Fig. ||). Which of two length scales 2d and I determines the total density of states will 
depend on the other parameters. 
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FIG. 2. The calculated distribution of the trajectory length in the F-layer: a) for different values of FS-transparency b) 
for small FS-transparency (T = 0.01). The double peak structure close to the smallest length originates from the first two 
reflections, whereas the distribution for long trajectories decays as exp(— Z/Z), with the average length I ~ 2d/T. 



Combining Eqs. (B|) and (jl4|) we obtain for the total sub-gap DOS 



N(E) = ^J2 E P(fc n )e 2marccosB/A , (16) 

cr— ± n— — oo 

where fc„ = 2n(E + ah) / 'vf ■ Similarly, from Eqs. (||) and ( E3 ) the total DOS for the energies above the gap is obtained: 

N(E) = I -£J2 E P(fc n )e- 2 l"l acosh l^ A l. (17) 

o"— ± n— — oo 

Thus, in both cases the density of states is fully expressed in terms of the Fourier transform of the trajectory 
length distribution. Most probably a real F-film has a non-uniform thickness due to the large scale roughness of the 
boundaries. Assuming a smoothly varying thickness we can take this into account by averaging expressions (ttq) and 
( |l7| ) over a Gaussian distribution of the thickness around a mean value d. This will also leads to a smoothening of 
the sharp features in DOS resulting from the lower cutoff in p(l). The qualitative behaviour will however not change. 
In our calculations, we have taken a width of the distribution to be of order 10% that correspond to condition of the 
experimentaL-3. 

IV. RESULTS AND DISCUSSIONS 

Eqs. (plf) and (O) express the DOS of a F-layer contacted by the superconductor in terms of the trajectory length 
distribution. Depending on the relative values of A, h, the Thouless energy vp/2d and T, the resulting DOS has 
different behaviors. We will concentrate mainly on the limit T<1. For the length dependence it useful to distinguish 
between normal metal with h = 0, a weak ferromagnet h ~ A of h < A and h > A and a strong ferromagnetic film of 
h ^> A . We analyze the DOS in each case for different values of drA/vp (2dr = I), being the relevant length scale 
in the limit T< 1. In the end, we summarize all results in a map of the DOS depending on drA/tip and dxh/vF- 

A. Normal film 

Let us start with a normal metal film (h = 0) contacted by the superconductor. The DOS is shown in Fig. [3|a-c, for 
different values of drA/v-p. In the limit of a very thin layer with dyA/up <C 1, the DOS has essentially the form of 
superconducting DOS with sharp peaks at E = ±A and zero DOS for energies inside the gap. By increasing dxA/vp 
the peaks are getting broader and a finite DOS appears at small energies. There is still an energy interval around 
E = with zero DOS (see Fig. |]a). Increasing dxA/T further leads to a suppression of the superconducting features 
of the DOS. The zero DOS interval become smaller and the DOS at other energies tends to be closer to the DOS of the 
normal state. Thus, the density of states develops a minigap around the Fermi level, which decreases with increasing 
dr roughly as vp/dr — vpT/dhi (£i m p/d) (see Fig. p3b). This minigap is related to the mean length of the trajectories 
I, which has a finite value, if d/£i mp is finite. The presence of weak bulk disorder in the normal film suppresses long 
trajectories. Formally, this was included in the distribution of the trajectory length as the exponentially decaying 
factor in Eq. (pd|), which leads to the finite mean length I = 2dlg£i mp /d. This act as an effective upper limit of order 
I of the length of the trajectories, which gives a lower bound to the energy of the, Andreev bound states. Similar 
features were found before within a tunneling modela and in the diffusive modclsHEHI in a disordered normal layer 
contacted by the superconductor. 

The peaks of E± = A originate from the first peak in the distribution p(l) at I ~ 2d, which at higher c?t A/vp move 
to lower energies given roughly by ±vp/2d (see Fig. pp). They originate from Andreev peaks (AP) resulting from the 
trajectories with I ~ 2d. By increasing dxA/vF the first AP moves to lower energies, and, when drA/wp ^ 1/T", the 
next AP appear at E± — A (Fig. [3p). In this case the DOS is close to the normal states values. Small deviations 
proportional to T display many AP, as shown in Fig. Be. Small peaks close the main peaks, which are more pronounced 
for the first AP, correspond to the second peak of p(t) at I ~ Ad. 






FIG. 3. The DOS vs. energy of the normal layer (h — 0) for different values of drA/vF ■ a) Suppression of the superconducting 
features of DOS with increasing drA/vF and appearance of a minigap of order vf/cIt- b) The first Andreev peaks (AP) at 
roughly ±«f/2g( correspond to the first peak in the distribution of the trajectory length, c) For very large cItA/vf the DOS 
has many AP, leading to small deviations of order T from the normal state DOS. 



B. Very weak ferromagnet, h < A 

Now we study the effect of the spin splitting in a ferromagnetic film on the DOS. First we consider the case of a 
weak ferromagnet, where the exchange field is of order of the superconducting gap, but h < A. A qualitative picture 
of the influence of h on the DOS follows from the condition m) for the formation of Andreev bound states in the 
ferromagnetic layer. This is the semiclassical quantization condition for coherent superposition of two subsequent 
Andreev reflection of a quasi-particle from S, which propagate along a trajectory of length I. It includes the phase 
gained by a quasiparticle of spin a along a trajectory (E + ah)l, and the phase shift produced by an Andreev reflection, 
arccos (E/A). If the energy of the quasi-particle is not close to A the phase shift resulting from Andreev reflection is 
different from zero. To obtain constructive interference the total phase must be an integer multiple of 27r, as follows 
from Eq. (g). The existence of an upper limit on the length of the trajectories (as discussed above), leads to the 
formation of a zero DOS interval (minigap) around E = ±h. Therefore, the total subgap DOS should be similar to 
the average of two by ±h shifted normal spectra. 

In column (I) of Fig. the DOS of the ferromagnetic film with h — 0.5 A is shown for different values of dxA/vp. 
Shifting the minigap leads to minima at E — ±h. The zero energy DOS becomes finite and increases with increasing 
drA/vp. At higher drA/vp ~ 1 the DOS has a smooth peak at zero energy and two dips at E ± h (see Fig. 0-Ia). 
Here, the AP are located at E ± A. By increasing dxA/vp the width of the dips decreases roughly as vp/dr and the 
first AP moves to energies below the gap, i. e., ~ ±vp/d (see Fig. |J-Ib). We can distinguish two domains of energies 
below the exchange field \E\ < h and energies above the exchange field \E\ > h. At higher drA/vp the first AP moves 
from \E\ > h to \E\ < h and the next AP appears at l-El > h. For the region \E\ < h the DOS shows a zero energy 
peak, if the first APs merge at E = (see Fig. S-Ic). This is a zero energy Andreev peak (ZEAP), which originates 
from phase shifting caused by the exchange field. Additional shifting of the AP results in an oscillatory behavior of 
the DOS in the domain \E\ < h. As shown in Fig. 0-Id, in the limiting case of drA/vp ~^> 1/T the DOS is close to 
the normal state value, exhibiting small deviations, which are of the order of the FS-interface transparency T. The 
deviations have the form of small oscillations at all energies. 



C. Weak ferromagnet, h > A 



In the case of ft, > A the suppression of the superconducting features of DOS occurs at lower dxA/vp, compared 
with the previous case of h < A (see Fig. ^-Ila). The subgap DOS has similar features as it had for energies below h 
in the case of h < A. The zero energy DOS increases to the normal state values as dxh/vp becomes of order unity. 
Then, the DOS has a smooth peak at E = and minima at E ± A, as is shown in Fig. Q-IIb-c. The AP at E ± A 
move to lower energies at a higher values of dxA/vp and form ZEAP, when they merge at E = (see Fig. H-IIc ). 
The size of the ZEAP is of order of T. 

The main feature of DOS in energies above the gap consists of sharp peaks at E — ±h (see Fig. H-IIb-d ). These 
peaks originate from a resonant transmission of the quasi-particles through the superconducting potential (A = 0) 
inside the F-film. A quasi-particle incident from the superconductor to the F-film with energy above the gap, is 
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FIG. 4. Density of states in a ferromagnetic layer in contact to a superconductor. Different columns labeled (I)-(III) 
correspond to the values of the exchange field given above. The thickness of the layer is increased from the bottom to the top. 
For explanation of the various regimes, see the text. 



scattered by the potential well whose width is determined by the length of the quasi-particle trajectory I inside F. 
The phase gained by the quasi-particle of spin a inside the potential well is given by (E + <rh)l. At E = —ah the 
incident and transmitted quasi-particles interfere constructively, which tesds to a reflection-less transmission. Similar 
effect were found before in normal metal-superconductor hetrostructuresE^I. For drh/vp > 1/T the subgap DOS shows 
an oscillatory behavior around the normal state value as a function of dh/vp (see Figs.H-IIc and 0-IId). The period 
of the oscillation is ir/2 and the amplitude is of order T. The amplitude is damped in the limit of drA/vp ^S> 1/T. 
Note that the phase of the oscillation depends on the energy. 

D. Ferromagnetic film 

Now let us consider rather strong exchange fields h^> A. In column III of Fig. [| we plotted N(E) for h = 10A at 
different dxA/vp. As shown in Fig. Q-IIIa the suppression of the superconducting features from DOS by increasing 
drA/vp is faster than in the weak ferromagnetic cases. In fact the DOS reaches almost to the normal state value for 
drh/v-p ~ 1. As long as drh/ vp^l, mainly long trajectories of I ~ I contribute to the zero energy DOS, which are 
well approximated by an exponential distribution of the form exp (—1/1). Replacing this approximation form of p(l) 
into the general expression of DOS Eq. (g) we find, that the zero energy DOS increases with drh/vp roughly as 

m = N w_z»tpM t (18) 

hi 1 — exp (—irvp/hl) 

This results is also applicable for the weak ferromagnetic case discussed above. As before, the shifting of the AP to 
lower energies (see Fig. 0-IIIb) leads to the formation of ZEAP at the Fermi level as is seen from Fig. 0-IIIc/d. Then 
the DOS develops coherent oscillation as a function of dh/vp with the period 7r/2. The amplitude and the sign of the 
oscillation depend on energy. Maximal amplitudes of opposite sign always occur at zero energy or at the gap energy 
(see Fig. [|-IIId).|-This results in an inverted energy dependence of the DOS by changing d, which has been observed 
in the experimentEJ. We have shown inE3 that our results is in a quantitative agreement with the experimental data. 

E. Maps of the proximity DOS 

Summarizing the above analysis we present a map showing the dependence of DOS on dxA/vp and dxh/vp for 
small FS-intcrface transparency. This map is shown in Fig. |5|. Various regions in the map are distinguished by 
different ranges of h/A and drA/vp (or equivalently dxh/vp). Along the diagonal lines h is equal to A and moving 
upwards dxA/vp (and consequently drh/vp) increases. The quarter circles are curves with constant dx, along which 
the ratio h/A is varying. In the following we discuss different regions according to this classification. The normal film 
corresponds to the vertical axis (h = 0), which consists of three parts. The first part is limited by drA/vp<l. Here 
the superconducting features are dominant at lower drA/vp and suppressed for drA/vp ~ I, showing a minigap at 
the Fermi level. The second part is limited by l>dxA/vp< 1/T, where the main feature is the first AP at energies 
~ ±vp/2d and a minigap of order vp/dx- Close to the boundary drA/vp ~ 1/T the second AP appears in the DOS. 
Finally, the third part is the region dxA/vp > 1/T, where the DOS contains many AP, appearing as small deviations 
(proportional to T) from the normal state DOS. 

In the case of a ferromagnetic film we distinguish the following regions: 

i) The strongly superconducting region is limited by the smallest quarter circle, in which the superconducting 
features are dominate the DOS. At non-zero h the zero energy DOS appears at larger drA/vp and the DOS 
increases to the value of the normal state in the domain close to the second quarter circle boundary. For the 
part h < A, there is a smooth maximum at E = between two minima at E = ±h. For h > A we have only a 
smooth peak at E = 0. 

ii) The intermediate regions limited by two quarter circles. In h < A part the main features in DOS is the existence 
of two dips at E ± h and the first AP. Close to the second boundary we observe a separation between two energy 
domains \E\ < h and \E\ > h. While the second AP peaks appears at energies above h, the first AP peaks move 
to energies below h domain. In the h > A part the DOS has a smooth peak at E = 0. In the domain close to 
the diagonal line h> A, there are also two resonance sharp peaks at E ± h which disappear in h 3> A regions. 
In both cases of h < A and h > A, shifting of the APs leads to the formation of a ZEAP. This happens at 
regions close to the third quarter circle. 
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FIG. 5. Map of the proximity DOS in a ferromagnetic layer showing the dependence on dr A/vy and drh/vp for a small 
FS-interface transparency T<1. In the case of high transparency T ~ 1 the region between 1 and 1/T is absent. 



iii) In the region above the third quarter circle the DOS is close to the normal state value (flat). There are, however, 
small deviations proportional to T, which have different origins in the different domains. In the limiting domains 
of ft ^C A and h ^> A (close to the respective axis') they consist of many AP and the oscillations, respectively. 
In the ft < A part there are many AP above ft and an oscillatory variation below ft. In the ft > A part and for 
weak exchange fields in subgap part we have oscillatory variations with an energy dependent sign and amplitude, 
which result from a collective shift of many AP. As the intermediate region the resonance peaks are present at 
E = ±ft. In the strong exchange field part these peaks are disappear. Here the oscillations of subgap DOS are 
produced by the first AP only. 

The DOS map for the case of high transparency T ~ 1 is similar to Fig. pi. The difference is the absence of the 
region between 1 and 1/T on both axis. In the remaining regions we have features similar to the T <C 1 case. The 
value of the minigap and the energy of first AP arc the same order. All features and variations including AP and 
oscillations are more pronounced than the T <C 1 case. 

V. CONCLUSION 

In summary, we have studied theoretically superconducting proximity effect in ballistic ferromagnetic layers. Within 
the quasi-classical formalism, we have obtained expressions for the density of states at all energies in a ferromagnetic 
metal in contact with a superconductor which are completely specified by Fourier transform of the length distribution 
of classical trajectories in the ferromagnet. The length distribution of trajectories depends only on the geometrical 
properties of the attached ferromagnet and the connected interfaces. Thus, the obtained expressions are applicable 
for ballistic FS structures of arbitrary geometry. We have calculated the length distribution for the film geometry in a 
quasi-ballistic model taking into account finite transparency of the FS interface, roughness at the film boundaries and 
weak bulk disorders. The density of states exhibits variety of structures depending on the values of the superconducting 
energy gap A, the ferromagnet exchange field ft, the thickness of F film d and the FS interface transparency T. We 
have observed many interesting features, like splitting of the subgap density of states for spin up and spin down 
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quasi-particles, a zero energy Andreev peak, resonant transmission peaks above the gap at E = ±h and oscillations of 
the DOS as a function of the film thickness. These effects have been explained in terms of the phase shift of Andreev 
bound states, caused by the exchange field. We have analyzed the density of states in the full parameter range and 
summarized the results in a map, shown in Fig. |5|. 
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